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Summary. — Sytematic uncertainties to the single top production cross section measurement 
at the ATLAS experiment has been studied. Different sources of systematic uncertainties such 
as detector luminosity, jet energy calibration, SM background normalization, PDF parameter- 
ization and others have been considered. Large scale Monte-Carlo events simulation has been 
performed to estimate the contribution of each source in the overall uncertainty. The study was 
done for cut based analysis as well as for multivariate analysis of the single top measurements 
in ATLAS. The total systematic uncertainties of the single top cross-section measurements in 
its three production channels have been estimated at 1 fb _1 and 10 fb _1 integrated luminosity 
offheLHC. 

PACS 14.65.Ha - Top quarks. 



1. - Introduction 

The single top production at the LHC in the three different channels (f , s and Wt) is expected 
to have a sizable production rate. This should allow to accurately measure these weak processes 
and determine the cross section in the different channels (TJ. As the detailed cut based and 
multivariate analysis show, the single top production processes in ATLAS could be discovered at 
the level of few hundreds of inverse pb (depending on the channel) of the integrated luminosity 
of LHC ID. 

However, several anticipated sources of systematic uncertainties will significantly affect the 
precision of this measurement. Detailed study of the systematic uncertainty sources has been per- 
formed within the ATLAS top physics working group. Influence of the uncertainties of the LHC 
luminosity, trigger efficiencies, jet energy calibration, b-jet tagging efficiencies, lepton identifi- 
cation, background normalization, initial and final state radiation (ISR/FSR), parton distribution 
functions (PDF) and b-fragmentation parameterizations on the single top cross section measure- 
ments were estimated using detailed Monte-Carlo studies of the processes in the ATLAS detector. 
Systematics were determined for cut based and multivariate analyses at different reference inte- 
grated luminosities. Estimation of the systematics from some uncertainty sources would have 
required very large statistics Monte-Carlo simulation. In those cases a re-weighting method was 
implemented, validated and used for estimation of the uncertainties. The re-weighting approach 
was necessary for PDF and b-fragmentation parametrization uncertainty studies. 
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2. - Parton distribution functions 

The PDF parametrization is based on a fit to a large number of experimental data points 
and inherits uncertainties from the data used in the fit. Then the PDF parametrization uncer- 
tainty can propagate to the uncertainty of the measured physical observable, for example, of the 
cross section. In order to estimate the resulting PDF uncertainties on the single top cross section 
measurements, first the PDF uncertainty on the selection efficiencies of the signal and main back- 
ground (f F) has been estimated. Then, the selection efficiencies were translated to uncertainties in 
the different single top channel cross section measurements. PDF uncertainties were estimated at 
the NLO level. For this, the MC@NLO Monte-Carlo generator together with the CTEQ6. 1 NLO 
PDF error set has been used for the event generation (for the Wt-channel AcerMC generator has 
been used instead). Events were generated at the so called central value PDF (corresponding to 
the PDF parameterization fit parameters at the global minimum of the fit) of the CTEQ6.1 |2| 
set. Each event was assigned a set of weights determined as follows, 



± ^ /i(*i,6;3f)-/2(* 2 ,g;Sf) 
Wi fi(xuQ;So)-f 2 (x 2 ,Q;So) 

Here, fa and fa are PDF values for a given hard scattering process (characterized by flavors 
/ and momentum fractions x of the initial partons, and by Q, the event scale) evaluated for the 
i th error PDF pair Sf . The terms in the denominator are obtained with the central value PDF. For 
the uncertainty calculation, events were taken into account together with the corresponding PDF 
weights. This means that weighted and normalized sums of selected events were determined as 
given by the formula: 



Lk=i 



where, M and m are the total and selected numbers of single top events respectively. The 
Wj^'s form a set of normalized 2N weight sums, which replace selection efficiencies e'"" x — > 
max(W ; + ,W ; ~) and e™" — > min(W ! + , Wj~) in the master equation for the up/down systematic 
uncertainties due to PDF, 



A+e = J£max(£ ( m£ "-£°,0) 2 and A"e = Jjjmx(e° - e™",0) 2 . 

In these equations, N is the number of global fit parameters and e° is evaluated at the central 
value PDF. Figure 1 shows the central value corresponding selection efficiency for the t-channel 
and for it (applying the t-channel selection) in the e - final states and the variation of the re- 
weighted selection efficiencies. 



3. - b-quark fragmentation 

The b-quark fragmentation determines the spectrum of the b-flavored hadron distribution as a 
parameterized function of a Lorenz-invariant energy-momentum variable - z. The b-fragmentation 
parametrization uncertainty affects the b-tagging efficiency and the b-hadron energy spectrum in 
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Fig. 1. - Variation of PDF re-weighted selection efficiencies. Left: t-channel into e -final state; right: tt 
(with t-channel selection) into e ± - final state 



Monte-Carlo simulation and therefore, it is an additional source of systematic uncertainty in 
the single top selection efficiency. The uncertainty was determined as a difference between the 
single top selection efficiencies due to the choice of the different b-fragmentation models, the 
Peterson [ 3 1 compared to the Lund-Bowler J5J (H model. The latter leads to an increase of the 
expected b-tagging efficiency (0.6) by 0.01 1. This result was used to estimate the corresponding 
selection efficiency uncertainty using an event weighting method. The normalized weight sum 
for selected events was determined considering the specific selection cut requirements on b-jets 
in the different single top channels: 



vf ± (s-channel) = -*^ 



1 " />?- 6±001 1 (exactly 2 truth b-jets out of M with p T > 30GeV are tagged) 
N r-j Pf 6 (exactly 2 truth b-jets out of M with p T > 30GeV are tagged) 



. 1 " P°' 6±o on (at least 1 truth b-jet out of M with« r >50GeV is tagged) 

w (t — channel) = — * ) — — jrr 

N *-[ ^ U b (at least 1 truth b-jet out of M with pr > 50GeV is tagged) 

where, Pf 6 and f> a6±0011 are the binomial probabilities for the z'-th event calculated at the 
b-tagging efficiency values: 0.6, 0.61 1 and 0.589. M is the number of true b-jets in the z'-th event 
satisfying a pr requirement. N and n are the number of total and selected events respectively. 
The accumulated normalized weight sums were used to calculate the corresponding up and down 
uncertainties on the selection efficiency. This is done using a similar procedure to the PDF 
uncertainty calculations. 

4. - Estimation of the total uncertainty 

The results of the PDF, b-fragmentation and all other systematic uncertainties were used 
to determine an overall uncertainty of the cross section measurement using the master formula, 
a = N '°^> ; Where, N tot is the sum of signal and background events provided by the Monte Carlo 
simulation, B is the sum of all background contributions, £ is the signal selection efficiency and 
££ is the luminosity. The effect of each source of uncertainty on these variables were combined 
and propagated to the measured cross section using a toy Monte Carlo method, which randomly 
generates Ntot according to a Poisson distribution, and randomly varies B and £ for every system- 
atic source by an amount chosen around its central value, according to a Gaussian distribution. 
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Fig. 2. - Single top cross-section measurement systematic relative uncertainties 



This procedure is performed a few thousand times and the RMS of the resulting distribution is 
interpreted as the total uncertainty. The total uncertainty is dominated by systematic effects due 
the background normalization. Reducing the background contamination to increase S/B would 
help to reduce the systematic uncertainty. Multivariate analyses were used to further optimize the 
selection beyond what can be achieved with cuts on the existing variables [ 1 1. Figure 2 shows 
the cross-section systematic uncertainties for the different channels at the 1 fb -1 and 10 fb -1 
integrated luminosities of the LHC calculated for the multivariate analyses. Table 1 provides the 
statistical and systematic uncertainties for these channels for the above integrated luminosities. 



channel 


Analysis at 1 fb 1 
± stat. % ± syst. % 


Analysis at 10 fb 1 
± stat. % ± syst. % 


t 


5.7 


21.7 


1.8 


9.8 


s 


64.0 


95.0 


20.0 


48.0 


Wt 


20.6 


48.0 


6.6 


19.4 



TABLE I. - Cross section measurement relative uncertainties for the different channels of the single top 
production at ATLAS. Multivariate analysis results 
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